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ABSTRACT 

We conduct two-dimensional axisymmetric (referred to as 2.5D) hydrodynam- 
ical numerical simulations of bubble evolution in clusters of galaxies. We inflate 
bubbles using slow, massive jets with a wide opening angle, and follow their evo- 
lution as they rise through the intra-cluster medium (ICM). We find that these 
jet-infiated bubbles are quite stable, and can reach large distances in the cluster 
while still maintaining their basic structure. The stability of the jet-infiated bub- 
ble comes mainly from the dense shell that forms around it during it's inflation 
stage, and from the outward momentum of the bubble and the shell. On the 
contrary, bubbles that are inserted by hand onto the grid and not inflated by 
a jet, i.e., an artificial bubble, lack these stabilizing factors, therefore, they are 
rapidly destroyed. The stability of the jet-inflated bubble removes the demand 
for stabilizing magnetic fields in the bubble. 

Subject headings: (galaxies:) cooling flows galaxies: clusters: general galaxies: jets 



1. INTRODUCTION 

Many clusters of galaxies harbor bubbles (cavities) devoid of X-ray emission, e.g., 
Perseus (Fabian et al. 2000) and Abell 2052 (Blanton et al. 2001). These low density 
bubbles are inflated by the jets launched by the active galactic nuclei (AGN) sitting at the 
centers of these clusters. 

The estimated ages of most bubbles (Birzan et al. 2004) is about one order of magnitude 
larger than the characteristic time of the Rayleigh- Taylor (RT) instability to destroy them. 
This observation has prompted many authors (e.g., Briiggen & Kaiser 2001; Kaiser et al. 
2005; Jones & De Young 2005) to invoke an ordered magnetic field at the edge of the bubble, 
or to consider the effects of viscosity (Reynolds et al. 2005), to stabilize the bubble against 
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the RT instability. Indeed, numerical simulations of non-magnetic bubble evolution show 
them to be disrupted quite rapidly (e.g., Briiggen 2003; Briiggen & Kaiser 2001; Jones & De 
Young 2005; Pavlovski et al. 2008; Reynolds et al. 2005; Robinson et al. 2004; Ruszkowski 
et al. 2004a, 2004b, 2008). Adding magnetic field make the bubble more stable (e.g., Jones 
& De Young 2005; Robinson, K. et al. 2004; Ruszkowski et al. 2007). However, it 
is not clear if magnetic fields can indeed supply the required stability (Ruszkowski et al. 
2007). In the simulations cited above the bubble where injected at off-center locations by a 
prescribed numerical procedure. This is of course not the way bubbles are formed in clusters. 
Evidence suggests that these bubble are formed by jets. We term these type of bubbles 
that were inserted numerically artificial bubbles. Different in that respect are the tower jet 
model simulations conducted by Nakamura et al. (2006; also Diehl et al. 2008), which did 
follow the evolution of a jet and a bubble. Nonetheless, we find some problems with this 
model, e.g., the angular momentum that is assumed in the jet is too large. 

Jets can inflate bubbles if their opening angle is large (i.e., wide jets), or if they are 
narrow but their axis changes its direction (Soker 2004, 2006; Sternberg et al. 2007; Stern- 
berg & Soker 2008). The change in direction can result from precession (Soker 2004, 2006; 
Sternberg & Soker 2008), random change (Heinz et al. 2006), or a relative motion between 
the ICM and the AGN (Loken et al. 1995; Soker & Bisker 2006; Rodriguez-Martinez et al. 

2006) . The outcome of a wide jet and a rapidly precessing jet is the same (Sternberg & 
Soker 2008) . For numerical reasons we conduct a study where we inflate bubbles with wide 
jets, and not with precessing jets. We use jets with Mach number of 10. Much faster jets 
form elongated bubbles, rather than 'fat' bubbles, i.e., bubbles with axes ratio close to unity. 
In addition, much faster jets cause back flow of hot gas that fills the region between the 
two bubbles, and therefore lead to the formation of one elongated bubble (Sternberg et al. 

2007) . In this Letter we limit our-self to comparing the evolution of jet-inflated bubbles (i.e., 
bubbles infiated solely by jets and not by a prescribed numerical procedure) with artificial 
bubbles, in the goal of emphasizing the necessity to inflate bubble self-consistently. For that, 
the description of the numerical code is brief, and we present only one case out of the several 
we simulated. 



2. NUMERICAL METHOD AND SETUP 

The simulations were performed using the Virginia Hydrodynamics- 1 code (VH-1; Blondin 
et al. 1990; Stevens et al. 1992), as described in Sternberg et al. (2007). Radiative cooling 
was not included. We study a three-dimensional axisymmetric flow with a 2D grid (referred 
to as 2.5D). We simulate half of the meridional plane using the two-dimensional version of 
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the code in spherical coordinates. The symmetry axis of all plots shown in this paper is 
along the x (horizontal) axis. Due to the fact that this is a Letter, in order to minimize the 
space required to show the plots, and due to the fact that the results of both halves of the 
simulated domain were either identical or exhibit the same large scale and stability behavior, 
we exhibit only a quarter of the meridional plane (half of the simulated domain) . 

In contrast to our previous papers, where gravity was omitted, we added gravity to the 
simulations. We assume a dark matter halo with a density profile as is given in Navarro et al. 
(1996). The dark matter is not affected by the evolution of the baryonic matter, therefore, 
the gravitational potential is constant. 



^NFw{r) = iTirlScPcritG 



^ _ \n{r/rs + 1) 



(1) 



where = — is a scale radius, Vy is the virial radius, c is the concentration factor, 5c — 

3 ^ 2 

^ )^^^ ^ c , and Pcrit — is the critical density of the universe at 2; = 0. 

We set the initial density profile of the gas to maintain hydrostatic equilibrium (e.g., 
Makino et al. 1998), i.e.. 



b 



Pgas{r) = Pgasfie ^yy + I j , (2) 

where h = AnrgScPcritpnT'pG/kTy. — '^'^^j^^J^'' is the virial temperature, and My is the 
virial mass. The unperturbed ICM temperature is 2.7 x 10^ K. We use a 256 x 256 grid, 
evenly spaced in the azimuthal direction. In the radial direction, the grid was partitioned 
using a geometric series with a common ratio of 1.0015 allowing for better resolution in the 
inner part of the simulated domain. More detail will be given in a forthcoming paper where 
an extended parameter space will be explored. 

We simulate the evolution of a bubble inflated by a jet, and compare it to the evolution 
of bubble introduced manually to the grid (an artiflcial bubble). The jet is injected at a 
radius of 0.1 kpc, with constant mass flux of Mj = 5Mq yr~^ (per one jet) and a constant 
radial velocity of Vj — 7750 km s~^, inside a half opening angle of a = 70°. The total power 
of the two jets is E2j — 2 x 10^^ erg s~^. The jet was active for a period of Atj = 10 Myr, 
from t — —10 Myr until t — 0, when the jet was completely shut-off. 

The spherical artiflcial bubbles was inserted in its full size at t = 0. with its volume 
and (constant) density about equal to that of the jet-inflated bubble a,t t — 0. The initial 
bubble's radius and density were set to Rabo — 4.5 pc, and pabo — 10~^^ g cm~^, respectively. 
We made several tests. We ran cases with initial artiflcial-bubble density three times larger, 
and three times smaller, than p^bo- We also ran cases with numerical grid with half and twice 
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t he resolution (i.e., 128x128 and 512x512 respectively) compared to our standard resolution. 
We were witness to only small differences between the runs with different resolutions and 
initial bubble density. Namely, our results and conclusions are not sensitive to the numerical 
details. 

With the parameters used here the bubble temperature is ~ 10 times higher than the 
ambient temperature. We can use a somewhat higher jet velocity that will result in a lower 
bubble density (Hinton et al. 2007) and in a higher bubble temperature (see Sternberg et 
al. 2007 and Sternberg & Soker 2007). This will be examined in a future paper. 

3. RESULTS 

In Figure [1] we show the density and velocity map of the jet-infiated bubble at t = 0. 
There are some flow structures which characterize the jet-infiated bubble, and that are not 
present in the case of the artificial bubble. As we shall see below, these are crucial for the 
future evolution of the bubble. (1) There is a dense shell around the bubble. (2) The bubble 
and the dense shell have radial momentum. In particular, the shell's front is relatively dense 
and has an outward velocity. (3) There is a circular flow, a vortex, around the lowest density 
region in the bubble. (4) Although the jet was shutdown, the jet material along the path 
form the center to the bubble still exists. In reality, the jet is gradually shut down (even if 
on a short time scale), which result in low density material behind the bubble, hence, this is 
a real feature. 

In Figure [2] we compare the evolution of the jet-infiated bubble with that of an artificial 
bubble. As is well documented in many non-magnetic simulations of artificial bubbles (e.g., 
Churazov et al. 2001; Robinson et al. 2004; Reynolds et al. 2005; Jones & De Young 
2005), as the bubble rises a vortex from below forms a flow along the axis that penetrates 
the bubble from below. This flow leads to the destruction of the bubble. In the case of an 
artificial bubble this occurs over a short distance, and the bubble clearly loses its shape. 

With the more realistic jet-infiated bubble the situation is very different: 

1. RT instability modes are seen on the front of the artificial bubble. The jet-inflated 
bubble does not suffer from such instabilities at early times. This is a result of the 
outflow velocity of the dense shell in front of the bubble. This interface is stable during 
the inflation phase (Soker et al. 2002; Pizzolato Sz Soker 2006), and the stability is 
maintained as long as the low density bubbles does not support the dense ICM against 
gravity. 
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Fig. 1. — Density and velocity map for the jet-inflated bubble at t = 0, the time the jet is shut 
off. The density scale is on the right in logarithmic scale and cgs units. The arrows represent 
the velocity of the flow: O.lCg < Vj < 0.5cs (shortest), 0.5cs < Vj < Cs, Cg < Vj < Sc^, and 
5cs < Vj < lOcg (longest in this case). Here Cg = 775 km s~^ is the sound speed of the 
undisturbed ICM. 
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2. The jet-inflated bubble reaches a much larger distance in the cluster before it starts los- 
ing its shape, as compared with the artificial bubble. This is a result of the momentum 
of the bubble and the dense shell around it. 

3. The jet material that lags behind the bubble, partially fills the region along the sym- 
metry axis. Our numerical grid forces the flow to be exactly axisymmetric. In a more 
realistic 3D flow, we expect this region to be spread around the symmetry axis, such 
that in projection the a bubble will still be observed as a more or less spherical bubble. 
In the case of an artificial bubble this region is filled with ICM dense gas, and in pro- 
jection the bubble will appear as a torus. In addition, in a future paper we intend to 
follow Reynolds et al. (2005) and increase the viscosity. We expect the higher viscosity 
to reduce the flow of dense matter along the symmetry axis. 

4. It seems that the vortex inside the jet-infiated bubble (Fig. [T]) stabilizes the sides of 
the bubbles as it rises. We see no sign of instabilities there, neither RT nor Kelvin- 
Helmholtz. 

5. There is a low density- high temperature (high entropy) gas lagging behind the bubble 
in a disrupted flow. This gas is mixed with the ICM and increase its entropy. What 
we find here is a relatively efficient way to heat the ICM. 

We would like to state that our results also support past conclusions that simulations 
of feedback in cooling flow clusters have to be more realistic and incorporate, among other 
things, a more realistic jet (Vernaleo & Reynolds, 2006). The inflation of the bubbles by 
the jet, and not by predescribed numerical recipe, might be one of the missing ingredients 
needed to make these simulations more real. 



4. DISCUSSION AND SUMMARY 

When the realistic jet-inflation process of bubbles is considered there are two stabilizing 
processes. (1) During the inflation phase the bubble is RT stable because the bubble-ICM 
interface is decelerating and expanding (Soker et al. 2002; Pizzolato & Soker 2006). (2) 
The outward momentum of the bubble and the dense shell around it implies that the low 
density bubble does not need to support a dense gas above it during the outward motion to 
large distances. This implies that the interface is RT stable. The result of these processes is 
that the bubbles can rise to a large distances from the cluster's center while still maintaining 
their general structure, without the need to invoke stabilizing magnetic field. Therefore, it 
is crucial that in the study of low-density bubbles in clusters of galaxies, the bubbles will be 



- 7- 



self-consistently inflated, rather than introduced artificially. It seems also that the vortex 
within the jet-inflated bubble suppresses instabilities on the sides of the rising bubble. 

Our results show that the inflating jet, even if completely shut down, result in a high 
entropy gas extending from the center to the bubble. This high-entropy gas mixes with the 
ICM. This mixing is a relatively efficient channel to heat the ICM. As seen at late times in 
Figure El the volume filled by the 'fractal' high entropy gas is non-negligible. 

The inflation of large, more or less spherical, bubbles close to the center of the cluster 
(termed 'fat' bubbles) requires the jet to be slow, Vj ~ 10"^ km s^^ ^ c, and the mass loss 
rate to be relatively large, ~ 1 — 5OM0 yr^^ (Sternberg et al. 2007; Sternberg & Soker 2008). 
In these papers we already mentioned AGN observations that support such a high mass loss 
rate, and listed some previous theoretical studies based on such an outflow. We here add that 
recent observations suggest that in many clusters, cooling flow does exist as the radiative 
cooling is not completely prevented (O'Dea et al. 2008). Namely, gas radiatively cools to 
low temperatures. Revaz et al. (2008) conducted numerical simulations and use them to 
suggest that cold blobs of gas participate in the feedback heating in cooling flow clusters, as 
suggested by Pizzolato & Soker (2005). Over all, it seems that high mass loss rate in jets 
might occur quite often in cooling flow clusters, as suggested by the moderate cooling flow 
model (Soker & Pizzolato 2005). Another issue is the opening angle of the jet. We note 
that most observed jets are narrow. However, as discussed in Sternberg & Soker (2007), it is 
possible that much of the energy and mass reside a wider and slower outflow, i.e., the narrow 
jet is the center of a wider jet. Our prediction of a wide, or rapidly precessing, jet will have 
to be tested in the future. 

Putting the results of this paper on a broader view, the usage of slow massive jets, 
that are either wide of rapidly precessesing, can account for some basic properties of cooling 
flow clusters: (1) Recycling of gas that cool from the ICM and flows toward the center. 
(2) Formation of fat bubbles close to the center. (3) Allowing the bubble to rise to large 
distances (the results of this paper). (4) Efficiently transfer energy form the central accreting 
black hole to the ICM. 

We thank John Blondin for his immense help with the numerical code. We thank Chris 
Reynolds (the referee) for helpful comments, and Edward Pope for his insightful remarks. 
This research was supported by the Asher Fund for Space Research at the Technion. 
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Fig. 2. — The evolution of the jet-inflated bubble (left column) and the artificial bubble 
(right-column), at three times as indicated. In the evolution of the jet-infiated bubble the 
jet was active from from t = —10 Myr until t = 0, when the jet was completely shut-off. 
The artificial bubble was introduces as a spherical bubble with a constant density at t = 0. 



